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SUMMARY 


Flutter characteristics for yaw angles between 15° and 90° were 
determined experimentally for two types of corrugation- stiffened panels: 
those with weak twisting stiffness and those with strong twisting stiffness. 

By mounting the panels on a remotely controlled turntable, good definition 
of the flutter boundaries was obtained by rotating the panels into and out of 
flutter. Flutter tests were conducted at M = 2 and M = 1.6 in the Langley 
Unitary Plan Wind Tunnel. Before flutter testing, vibration tests and analyses 
were also performed. The experimental flutter data is compared with flutter 
theory for orthotropic panels utilizing quasi-steady aerodynamics. In total, 
five different corrugated panels were tested consisting of one single skin 
panel having a length- to -width ratio of 5 on clamped supports and four 
different square double skin panels on discrete flexible supports. The 
investigation indicated that flutter speed for corrugated panels is highly 
dependent on yaw angle. Reasonable flutter correlation between analysis and 
test was obtained for moderate yaw angles but extreme sensitivity to 
structural parameters made the correlation at large yaw angles uncertain. 



INTRODUCTION 


A candidate for the thermal protection system of high speed reentry and 
hypersonic vehicles is the corrugated metallic panel constructed of high 
temperature alloy. Such a panel may he idealized as an orthotropic plate with 
either flexible or rigid supports. Flexible supports are often necessary to 
accommodate thermal expansion without creating large thermal stresses. 

Ideally, such a panel would be aligned on the vehicle with the corrugations 
parallel to the air flow. However in a typical flight trajectory various 
degrees of cross flow will be encountered which, analysis indicates, may 
reduce flutter margins drastically (ref. l) . Since available wind tunnel data 
were insufficient to substantiate these analytical predictions, a test program 
was conducted to provide such data. Specifically, various types of orthotropic 
panels were designed, analyzed, fabricated, and tested in the Unitary Plan Wind 
Tunnel at the NASA Langley Research Center. The results of this test program 
and the comparison with theory are the basis of this report. 

Many individuals at Grumman contributed to the work reported herein. 

The author wishes to acknowledge the efforts of Mr. John Valentine in designing 
the panel models, Mr. Edward Leszak for supervising the manufacture of the 
panels, and Messrs. Edward Ham, Timothy Foley, Anthony Longano and Paul Chase 
for testing the panels. The author is also grateful to Mr. John Smedfjeld 
for advice and suggestions. 


NOMENCLATURE 


a panel length (x-direction, fig. 8) 

b panel width (y-direction, fig. 8) 
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panel bending stiffnesses in x- and y- directions, respectively 
panel twisting stiffness 
frequency 
flutter frequency 

deflectional, rotational, and torsional 
spring constants, respectively, per unit length 
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non-dimensional deflectional spring constant 
non-dimensional rotational spring constant 
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non-dimensional torsional spring constant 
Mach number 

number of half-waves in x,y directions, respectively 
pressure 

pressure of center of turntable 

dynamic pressure of airstream 

dynamic pressure of airstream at flutter 

Cartesian coordianates of panel (fig. 8) 


Vm 2 -1 

yaw angle 

dynamic pressure parameter, 
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dynamic pressure parameter at flutter 
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Poisson's ratio associated with curvature in y- and x- 
directions, respectively 


TEST APPARATUS 


Panels 

The test panels (aluminum) were designed to simulate various stiffness 
parameters and support conditions. The summary of panel types is given in 
table 1. The single skin panel is representative of the torsionally weak 
construction, and the double skin panel of the torsionally strong construction. 
The single skin panel was formed from a 0.0254 cm (0.01 in.) sheet, and was 
bolted directly to a 1.27 cm (0.5 in.) thick aluminum supporting plate as 
shown in figure 1. The bolt spacing was 2.54 cm (l.O in.). The double skin 
panels were formed from 0.0203 cm (0.008 in.) sheets, and joined by spot 
welding along their lengths at the flat part between corrugations. Except for 
the different cross-sections, the double skin panels were identical in 
construction and were similarly supported. These panels were bolted directly 
to the circular supporting plate on two opposite sides and mounted on discrete 
flexible clips on the other two sides, as shown in figure 2. The clips were 
attached to the supporting plate (fig. 3) with either both the inner and outer 
screws or with just the outer screws. This provided two different support 
stiffnesses. Since the clips were manufactured in two thicknesses, 0.0595 cm 
(0.024 in.) and 0.0457 cm (0.018 in.), a total of four different spring 
stiffnesses were available. In order to evaluate the effect of beading on 
panel flutter the double skin panel with the flat cover sheet was tested twice, 
once with the flat side exposed to the flow and the second time with the 
corrugations exposed to the flow. 

As shown previously, the test panels were mounted on a circular plate 
which was attached to a turntable. For the double skin panels, fairings were 


4 



provided to smooth the air flow over the panel. The single skin panel 
corrugations were closed out at the ends and thus this panel did not require 
fairings. The turntable itself was mounted in the splitter plate which was 
projected into the airstream from the tunnel sidewall so that the tests were 
conducted free of the tunnel boundary layer. The whole arrangement (with a 
double skin panel installed) is shown in figure 4. The cavity pressure behind 
the panel was controlled manually by a 2.54 cm (l.O in.) diameter line attached 
to a vacuum pump. 


Instrumentation and Data Acquisition 

Each panel specimen was instrumented with eight single-arm strain gages 
and five iron -cons tan tan thermocouples. The locations of the strain gages 
and thermocouples are shown in figure 5 for the single skin panel and for a 
typical double skin panel. In addition, two flexible supports of each 
thickness were instrumented with strain gages. 

Signals from the strain gages were used during testing to detect the onset 
of panel flutter and to measure flutter frequency. They were monitored on an 
oscillograph and a dual beam oscilloscope. The signals from the thermocouples 
were recorded on a Brown multi-point recorder, and observed during test to 
assure that temperature differential was a minimum during measurement of the 
flutter threshold. All monitoring and recording was done at the side of the 
wind-tunnel test section containing the viewing window, so that yaw angle could 
be measured visually and the observation of the specimen could be maintained 
throughout the tests. High speed 1 6 mm movies were also taken to record panel 
motion. 

In addition to panel instrumentation, a calibration plate (similar to 
panel supporting plates) with eleven pressure taps was installed in the 
turntable to measure pressure variation over the panel. Seven pressure taps 
were located along the centerline of the turntable in the direction of the 
flow. Spacing of these taps was 8.89 cm (3*5 in.). Pressures were recorded 
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on a scanivalve system for M = 2 and M = 1.6 and for the wind-tunnel range of 
dynamic pressures. 


Wind Tunnel 

Flutter tests were conducted in the Iangley Unitary Plan Wind Tunnel at 

Mach numbers 2 and 1.6. The dynamic pressure in that tunnel is continuously 

variable. The maximum levels attained at Mach numbers 2 and 1.6 were 85. 1 
P P 

kR/m (1779 psf) and 70.6 kR/m (1V75 psf), respectively. The wind tunnel 
static temperature varied between 3H Q K and 328°K (lOO to 130°F). 


TEST PROCEDURE 


Vibration and Static Tests 

Vibration surveys and static tests were performed to check panel and 
support stiffnesses. For the vibration surveys, the panels were mounted in 
the circular supporting plate and excited by one or two air shakers. Rode 
lines were determined using a deflectometer mounted on a movable arm. These 
tests were conducted without simulating the effect of the cavity. Once the 
panels were in the turntable on the splitter plate, the shake tests were 
repeated at ambient pressure to determine the effects of the actual cavity. 
Ro appreciable change in frequencies was observed. Clip flexibilities were 
also checked statically by loading the clip with forces and moments and 
measuring the resulting displacements and rotations. 


Flutter Tests 


Before the test models were inserted into the airstream, 
plate with the pressure taps was installed in the turntable. 


the calibration 
Steady-state 
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pressure distributions were measured over the turntable area where panel models 
were to be located. These pressures were used to determine pressure loading on 
the panels, and the cavity pressure setting which would result in zero average 
differential loading. The measured pressure distributions for four different 
dynamic pressures levels at M = 2 and M = 1.6 are shown in figure 6. All 
values are relative to the pressure at the center of the turntable. Low 
pressure variations over the panel were obtained at M = 2 at low and inter- 
mediate dynamic pressures. However, large pressure variations existed at 
M = 1.6 and M = 2 at high dynamic pressure levels. 

For the first three panels tested, flutter boundaries were obtained for 
Mach numbers 2 and 1.6. However, it was determined from the data that the 
effect of Mach number on flutter speed corresponded to the well known factor 
of l/P. Thus, tests at M = 1.6 were abandoned, since at M = 2.0 the tunnel 
was smoother and the cavity pressure could be better controlled. The usual 
test procedure consisted of establishing constant temperature, dynamic pressure, 
and Mach number, and then rotating the panel away from the stiffest direction 
until flutter was initiated. During this time, cavity pressure was manipulated 
to maintain zero average pressure loading on the panel. At a particular 
dynamic pressure, a flutter point was defined as the smallest angle at which 
flutter could be induced by changing cavity pressure. Once the flutter point 
was reached, the panel was rotated back toward the stiffest direction, and the 
dynamic pressure increased by a certain increment. This procedure was 
continued until the maximum wind-tunnel dynamic pressure was attained. In 
this way the flutter boundary as a function of yaw angle was obtained. 

Additional check points were obtained by lowering the dynamic pressure by 
increments and repeating the above procedure. Since the turntable was power- 
driven, quick excursions into and out of flutter were possible. A typical 
trace of a strain gage response before and after initiation of flutter is 
shown in figure J. 
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RESULTS AND DISCUSSION 


Experimental and Analytical Frequency Comparisons 

The combinations of panels and support stiffnesses studied are listed in 
table 2, and the corresponding measured and calculated frequencies are given 
in table 3. The coordinate system and the various panel and support parameters 
are defined in figure 8. Panel stiffness properties were computed using the 
method of reference 2, and clip stiffness properties were calculated using 
strain energy methods. Natural frequencies were then obtained using a Galerkin- 
type solution (ref. 3)» 

For panel 1, measured frequencies in the long direction, m, fall between 
the simple -clamped and clamped- clamped calculated values. This indicates that 
although the panel is essentially clamped along its long sides (x-edges), its 
short sides (y-edges) have a finite rotational restraint. Unfortunately, at 
the present time, the theory is incapable of analyzing flexible-clamped plates, 
and the actual degree of fixity was not determined analytically. In the short 
direction, the poor correlation is due to the difficulty in obtaining test 
modes and to the additional stiffness provided by the beads. 

For panels 2 through 5, except for the fundamental mode, the initial 
comparison of frequencies between test and analysis was poor. Since the 
measured individual clip deflection spring constant compared very well with 
the calculated value, the more doubtful values of rotational and torsional 
spring constants were changed to improve correlation. Clip stiffness properties 
corresponding to the "best fit" with the measured frequencies are listed in 
the last column of table 2, and the resulting frequencies in the last four 
columns of table 3- 

The improvement in frequency correlation is shown graphically in figure 9, 
where the calculated frequencies for the initial and "best fit" support 
stiffnesses and measured frequencies are plotted against weak direction mode 
numbers. The torsional spring constant lOp was found to be the only effective 



parameter influencing higher modes. The rotational spring constant 
influenced only low modal frequencies. 

As seen in figure 9(c), the initial calculated frequencies for panel 4 
(clips attached to the panel only by the two outer screws) compared reasonably 
to measured frequencies. This reinforced the assumption that the calculated 
panel stiffnesses were correct , and that for clips attached by three screws , 
frequency discrepencies were due to the underestimation of support stiffnesses. 

Panels 2 and 3 were identical except that different sides' were exposed 
to the flow. One would expect identical frequency responses. However, the 
vibration tests have indicated consistently different behavior, especially 
for second and higher modes and is reflected in the different "best fit" 
stiffness values of supports. These differences are not fully understood but 
are attributed to panel warpage and to slight variations in support attachments. 

Similar difficulty in natural frequency correlation for corrugated panels 
was experienced previously by other investigators (ref. 4). 

Experimental Flutter Results 

Flutter boundaries as a function of yaw angle for all panels listed in 
table 2 are presented in figures 10 through 15. Because of wind tunnel dynamic 
pressure limitations, no flutter points were reached between 0° and 15°, and 
in most cases flutter was obtained only between 30° and 90° • An unusually 
small amount of data scatter was encountered. This is attributed to: (l) the 
ease of initiating and stopping flutter by rotating the turntable; and (2) the 
small temperature effects due to the type of panel design, in which thermal 
expansion was accommodated by beads and discrete flexible supports. 

As can be seen in figures 10 through 15, yaw angle has a considerable 
effect on the flutter speed of each of the panels tested. For example, in 
figure 13, the rotation of panel 3 from 90° (flow in the weak direction) to 
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35° results in a threefold increase in flutter X. If the yaw angle is 
decreased further, greater increase in flutter speed may be realized. 

In addition to the panels listed in table 2, an attempt was made to 
obtain flutter boundaries for two isotropic panels having a length-to-width 
ratio of five. These panels, of 0.0305 cm (0.012 in.) and 0.0406 cm 
(0.016 in.) thicknesses, were mounted in the same supporting plate as used 
for the single corrugation panel. Wo consistent flutter points were obtained, 
however, due to the constantly changing temperatures which induced buckling 
of the panels. 

Other difficulties encountered during testing were with the cavity 
pressure control and pressure variations over the panels. At M = 1.6, and 
at M = 2.0 at low dynamic pressures, leakage around the turntable was too 
high to permit satisfactory manual control of the cavity pressure. For this 
reason, the normalized flutter dynamic pressures (q/(3) shown in figure 10 
are higher at M = 1.6 than at M = 2.0, despite the fact that reference 5 
suggests that the M = 1.6 boundary should be below the M = 2.0 curve. Also, 
due to higher pressure variations over the panels at high dynamic pressures 
(see figure 6), experimental flutter trends at low yaw angles were more 
difficult to define. This was especially true for panel 4, which had the 
lowest support stiffness. (See figure 1*0. 

Considering all of the above factors, the M = 2 results were judged 
superior to the M = 1.6 results, and accordingly the former were used for 
correlation with theory. 


Correlation with Theory 

The flutter analysis used for comparison with the experimental data is a 
generalization to account for flow angularity of the procedure due to 
reference 3. This is a Galerkin-type solution for an orthotropic panel on 
continuous deflect ional, rotational, and torsional springs along two opposite 
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edges, and simply supported along the other two sides. Pressure loading is 
obtained using static aerodynamics. Flutter is defined as the lowest dynamic 
pressure corresponding to a coalescence of any two modes. Numerical results 
presented here represent converged values using as many as thirty modes. 

At this point, it should be emphasized that the above analysis requires 
two opposite sides of the panel to be simply supported, whereas, for the test 
panels these sides were clamped. For the square panels (panels 2 through 5), 
this approximation is certainly valid, since the effective length of the panel 
in the y direction is extremely large, and the x edges do not influence either 
the panel natural frequencies or panel flutter dynamic pressures. For panel 1, 
which has four clamped sides, it was found that the best flutter results, at 
least for the range of yaw angle where test data was available, were obtained 
by making the short rather than the long sides simply supported. These results 
are presented in figure 11, where in addition, flutter points at A = 0° and 
90 Q for a clamped- clamped analysis (reference 6) are shown . At A = 90°, A 
is the same whether the y-edges are simply supported or clamped. 

For panel 4, for which initially good frequency correlation was obtained, 
good flutter comparisons between test and analysis are also evident (figure l4) 
for the very limited range of yaw angle. The deviation of test values at 
lower angles can be attributed, as stated previously, to the buildup of 
pressure loading at these dynamic pressure levels (see figure 6(a)). 

For all other panels, correlation is not as good. The discrete flexible 
supports introduce an additional set of variables having a large influence on 
flutter. Calculated flutter speeds at high yaw angles are particularly 
sensitive to the torsional stiffness of the supports. This is illustrated in 
figure 13, where the calculated flutter parameter A is plotted for K^, = 0.016 
and K^, = 0.30. = 0.016 is the value calculated using strain energy methods 

for one clip and dividing the value by clip spacing. K^, = 0.30 represents the 

"best fit" value based on test natural frequencies. The A at A = 90° is 

_ _ cr 

six times greater for = 0.30 than for = 0.016. Because of this wide 
range of possible calculated values, it was felt that a more realistic 
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approach to correlation would be to present both of the calculated flutter 
boundaries with the test results. Thus for panels 2, 3 and 5, two analytical 
flutter boundaries and one experimental boundary are given in figures 12, 13, 
and 15 , respectively. 

Test results for panels 3 and 5, in figures 13 and 15, fall between the 
two analytical flutter points at A = 90° and approach the upper boundary at 
lower yaw angles. Thus correlation improves at small and moderate yaw angles 
where the effect of torsional spring supports is smaller. Repeatability of 
test results is very evident for these two similar panels. It should be noted 
that for both panels corrugations are exposed to the flow. 

The configuration of panel 2 is identical to panel 3 except that the flat 
sheet is exposed to the air flow. The test flutter boundaries (figs. 12 and 
13) are, however, quite different. X values for panel 2 are less than a third 
of those for panel 3. Some of the difference can be explained by the different 
"best fit" stiffness characteristics as shown in table 2. However, it should 
be noted that panel 2 is the only one for which the test flutter boundary falls 
below the lower analytical boundary. This implies that panel beading is 
beneficial for panel flutter. 

A possible reason for the "best fit" flutter boundaries falling above the 
test data for large yaw angles is the overestimation of D for the double 
skin panels. Reference 2 indicates that the correlation of the calculated I ) 
with test values is highly dependent on the type of cross-section. Furthermore, 
cross-sections of the test panels in this investigation are open at the ends 
and accordingly their full torsional stiffness is not developed until some 
distance from the supports. 

Flutter frequencies were measured concurrently with the flutter dynamic 

pressures. The comparison of experimental and theoretical flutter frequencies 

is shown in figure 16 for panels 1, 2, k , and 5* As expected, good correlation 

was obtained for panel 1. For the other panels good correlation only exists 

where the test X is close to the calculated X . This reinforces the 
cr cr 
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contention that mode coalescence is the flutter mechanism. 


It should be noted that the analytical flutter point at A = 0° is the. 
coalescence of the first two modes in the strong (x) direction, and as such 
represents perfect alignment of the air flow with the corrugations. Any 
deviation will involve modes in the weak (y) direction ( m = 1, n = 1,2,3), 
and will lower the flutter speed and flutter frequency. Therefore, it is 

unrealistic to expect to reach experimentally the analytical X at A = 0° 

D cr 

for orthotropic panels with large 12 ratios. 

D„ 


CONCLUSIONS 


Flutter characteristics of two types of corrugation-stiffened panels at 
various yaw angles were studied both experimentally and analytically. Flutter 
tests were made at Mach number 2 and 1.6. Analyses were based on quasi-steady 
aerodynamics. In addition, vibration tests and analyses were performed. 

The experimental part of the study has shown that reliable panel flutter 
results can be obtained if care is taken in eliminating thermal stresses in 
panel models. A unique feature of this study was the use of a remotely 
controlled turntable to rotate test panels. Rotating the test panel to 
initiate and stop flutter clearly delineated the onset of flutter. The results 
indicate that flutter speed for corrugated panels is highly dependent on yaw 
angle . 

For moderate yaw angles, when a large component of the flow is along the 
corrugations, correlation with theory was good, but as the flow became 
perpendicular to the corrugations, large discrepancies between test and theory 
are apparent. Since for this configuration flutter is extremely sensitive to 
boundary conditions, a better definition of the structure is needed. 
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It can also be concluded that beaded surfaces have a stabilizing effect 
on flutter, though their exact influence is difficult to determine because 
of the large effect of other parameters. 
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Table 2 Panel Flutter Models 
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^Except for panel 1, x-edges are assumed to be simply supported. 




Table 3 Panel Natural Frequencies 



S-C = Simply supported - Clamped 
C-C = Clamped - Clamped 
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Figure 4 Double Skin Panel Mounted 
in the Splitter Plate 
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Figure 7 Strain Gage Response During Flutter 
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Figure 8 Schematic of Orthotropic Panel and 
Support System 


25 



• Measured 
A Calculated 
□ Calculated "Best Fit" 



12 341234 

Mode Number, n Mode Number, n 

(a) Panel 2 (b) Panel 3 



123412 34 

Mode Number, n Mode Number, n 

(c) Panel 4 (d) Panel 5 


Figure 9 Comparison of Experimental and Theoretical Weak-Direction Panel 
Natural Frequencies 
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Figure 11 Comparison of Experimental and Theoretical 
Flutter Boundaries for Panel 1 at M = 2. 0 




Figure 12 Comparison of Experimental and Theoretical 
Flutter Boundaries for Panel 2 at M = 2. 0 
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Figure 16 Comparison of Experimental and Theoretical Flutter Frequencies at M = 2 
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